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Giant enhancement of the nonhnear response of a single crystal Nb sample, placed in a pumping 
ac magnetic field, has been observed experimentally. The experimentally observed amplitude of the 
output signal is about three orders of magnitude higher than that seen without parametric pumping. 
The theoretical analysis based on the extended double well potential model provides a qualitative 
explanation of the experimental results as well as new predictions of two bifurcations for specific 
values of the pumping signal. 

PACS numbers: 74.25.Nf, 74.60.Ec 



I. INTRODUCTION 

It was well established in the sixties that nucleation 
of the superconductivity in magnetic field parallel to the 
sample surface begins when magnetic field H becomes 
lower than surface nucleation field i?c3 (-^cs ~ l-69-ffc2) 
The experiments at low frequencies reveal a non- 
linear nature of the ac response of the superconductors 
in surface superconducting state Q, |^ . A critical state 
model for the superconducting surface sheath was used 
for theoretical explanation of the experimental findings 
nil. Recently our study of the nonlinear properties 
of single crystal Nb showed that the critical state model 
hardly fits the new experimental data and a simplified 
two level model was used for the theoretical explanation 
of the experiment j5|. 

The above mentioned experiments in superconducting 
materials related to the nonlinear effects such as har- 
monic generation or rectification. But there is another 
class of nonlinear effects - parametric phenomena. Para- 
metric phenomena are important in fundamental physics 
and for applications. Parametric amplification of the 
weak microwave signals in Josephson junctions were stud- 
ied in detail in thew seventies (see and references 
therein). Experimental observations of parametric am- 
plification of the weak microwave signals were possible 
because of the nonlinearity of the Josephson junction Q. 
To the best of our knowledge parametric phenomena in 
surface superconducting states of bulk superconductors 
has yet to be studied. 

The present study is a continuation of our previous ex- 
perimental work and the theoretical analysis of the non- 
linear dynamics in the surface superconducting states of 
a single crystal Nb In these experiments the dc and 
low-amplitude high-frequency ac fields have been applied 
along the long axis of a sample. It was found that in the 
steady-state regime the nonlinear effects exist only in the 
surface superconducting phase. A unique aspect of this 
system is that the non-linearity does not exist a priory, 
but rather appears only in the presence of an imposed dc 
field. 

The new experiments have been performed in the same 
system, subjected to an additional low-frequency excita- 



tion field, and the amplitude of the output signal was 
measured as a function of the amplitudes of the external 
fields including that of the ac excitation field. In this 
experiment the nonlinearity of the system plays a double 
role. Usual parametric experiments deal with two signals. 
The system is excited by weak and strong (pumping) 
signals. The frequency of the pumping signal is twice as 
large as the frequency of the weak signal. Due to the non- 
linearity of the system, power from the pumping signal 
transfers to the weak signal 6j. In our experiment, the 
system was excited by an amplitude modulated ac field 
having some carrier and modulation frequencies. The 
spectrum of the ac field does not contain the modula- 
tion frequency. The modulation frequency in a system 
appears when the sample is in a nonlinear state. Only 
for a nonlinear system a pumping signal at a frequency 
twice as large as the modulation frequency affects the 
nonlinear response. The theoretical analysis, based on 
the extended double well potential model of two surface 
states, provides a qualitative explanation of the exper- 
imental data and predicts some new results which are 
partially confirmed by experiment. 



II. SAMPLE PREPARATION AND 
EXPERIMENTAL TECHNIQUE 

A rectangular (10 x 3 x 1 mm^) sample was cut by an 
electric spark from a single crystal niobium bar which had 
been fabricated by electron beam melting of high purity 
Nb precursor. The [100] crystalline direction is parallel 
to the long axis of the sample and the [010] direction is 
perpendicular to the sample plane. After mechanical and 
chemical polishing the crystal was annealed for one hour 
at T=1800 K under a vacuum of IQ-"^ Torr. These pro- 
cedures provided a high quality sample, as confirmed by 
the high resistance ratio Rsqok / Riqk ~ 300 and unper- 
turbed critical temperature Tc — 9.15 K. The details of 
the employed Nb crystal characterization, including H-T 
phase diagram, have been published elsewhere 0. 

The nonlinear response experiment was carried in a 
modified commercial SQUID magnetometer, also used 
for measurements of the dc magnetic properties. The Nb 
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FIG. 1: (Color online) Block diagram of the experimental 
setup. HFC - high frequency generator, LFG - low frequency 
generator, HPF and BPF - high and band pass filters, LCK - 
lock-in amplifier. 
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FIG. 2: (Color online) The amplitude of the output 
signalAn (H) as a function of the dc field H at different pump- 
ing amplitudes h2si. 



rod was placed in dc and ac magnetic fields applied along 
the [100] direction. The excitation ac field generated by 
a high frequency generator operating in a constant cur- 
rent regime had the form of an amplitude modulated 
ac field h{t) = ho{l + acos^t) coauit. Typical values 
of above parameters were < ho < 0.4 Oe, a « 0.9, 
uj/2n = 3.2 MHz, and n/2Tr = 1455 Hz. A parametric 
pumping field hp(t) = h2Q. cos{2VLt) was achieved using a 
copper solenoid of the commercial SQUID magnetometer 
system. 

A block diagram of the experimental setup is shown in 
Fig. n A 65-turn primary copper coil driven by a high 
frequency generator (HP3325A) produces an amplitude 
modulated ac field h{t). Signals with frequencies and 
20, were generated by an Agilent 33250A model gener- 
ator. The 1500 turn secondary copper coil was wound 
directly on the primary coil. The length of both coils 
was 15 mm. The voltage drop over the secondary coil was 
measured by a lock-in amplifier (EG&G 7265). The mod- 
ulation frequency $7 was used as a reference frequency for 
the lock-in amplifier. The ac field amplitudes /iq a-nd /120 
were measured by an additional small probe coil wound 
under the primary coil. For the sake of clarity a probe 
coil is not shown in Fig. ^ 

One should note that the ratio between the pumping 
frequency and the modulation frequency was equal ex- 
actly to 2 in our experimental setup. 

When the ac excitation is applied to a superconduct- 
ing sample in the nonlinear state, the magnetic moment 
of the sample oscillates at the harmonics of the funda- 
mental frequency lj, at frequencies ± fJ, and at the 
frequency of modulation SI and its harmonics. The sec- 
ondary coil converts these magnetic moment oscillations 
into an ac voltage signal. In the experiments we have 
measured, by means of lock-in detection, the amplitude 
of the signal at frequency as a function temper- 



ature, dc magnetic field 77, ac field amplitude and 
amplitude of parametric pumping /i2o . 



III. EXPERIMENTAL RESULTS. 

Here we restrict our consideration to the experimental 
data which allow a direct theoretical explanation. Fig. |21 
shows the field dependence of the rectified signal Aq_ mea- 
sured at r = 7.5 K for different amplitudes of /120 com- 
pared with the data obtained in the absence of a pumping 
field, /i2n = 0. 

As evident from Fig. |3 the application of parametric 
pumping leads to a giant enhancement of the rectified sig- 
nal. Under our experimental conditions the enhancement 
of reaches two hundred times, the signal in absence of 
parametric pumping. The present result confirms our re- 
sent observation Q that under stationary conditions the 
pronounced rectified signal appears only in the surface 
superconducting state. 

The measurements of An as a function of /120 for con- 
stant H and provides interesting information about 
nonlinearity of the sample. In figures |31 and ^ we present 
the amplitude of the output signal A^ as a function of 
the amplitude of the pumping field /i2o for different mag- 
nitudes of H and amplitudes of the dc and ac fields. 
Fig. 13 shows that increasing the amplitude of excitation, 
/iQ, does not always lead to a straight enhancement of the 
amplitude of the nonlinear response. Nonlinear systems 
with strong nonlinearity frequently demonstrate this type 
of behavior. 

The experimental data presented in Figs. 13 El are 
very complicated. In order to straighten out this many- 
parameter dependence a theoretical basis is needed which 
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FIG. 3: (Color online) The amplitude of output signal An as 
a function of pumping amplitude ^20- 
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FIG. 5: (a) Schematic representation of the free energy of a 
system as a function of the magnetic moment ; (b) a tilted 
double well system ; (c) a two-level system. 
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of them corresponds to a local minima in the free energy. 
For the theoretical description of the nonlinear response 
we use the well-known Landau-Khalatnikov equation of 
motion for the magnetic moment M of the surface state 



dM 



(1) 



where F is a phenomcnological coefficient which hereafter 
is assumed to be equal to unity, (F ~ 1), and the total 
Ginzburg-Landau free energy G is composed of a series 
of parabohc branches, G = ^G„, corresponding to the 
number n of the single flux quantum $0 iii the state Gn 



G„ = ^ (M - M„)2 
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FIG. 4: (Color online) The amplitude of output signal An as 
a function of pumping amplitude /i2n at different magnitude 
of magnetic field H . 



we will discuss in the next section. 



IV. THEORETICAL ANALYSIS AND 
DISCUSSIONS 

In the experiments described above the magnetic field 
H is the large constant field, and h [t) is a low amplitude 
ac field. It turns out that for a given value of H there 
are many superconducting states which can be labelled 
by the number n of fluxons pinned in this level, and each 



where M„ = ti^o/ttR'^, and e is the energy to nucleate 
a fluxon. 

According to Eq. ^ , the free energy as a function of an 
external field is composed of a series of parabolic branches 
(corresponding to different n) with centers located on an 
ascending curve (Fig. EJ. 

Instead of considering the free energy Eq. Q we re- 
strict ourselves to the relevant initial and final states. In 
line with this assumption let us replace the many-wells 
system shown in Fig. by a tilted double-well system 
(Fig. Eb), where the right minimum at M = corre- 
sponds to the surface superconducting nucleation field 
Hc3, and the left minimum is related to some magnetic 
moment M induced by an external field H. The distance 
between the minima is proportional to — H, 

Hence, the adopted form of the free energy 

AG AP p^M* 
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contains two parameters 

p^{Hc,-H)-' : 7 = ^ (4) 

J^C3 

Note that in our previous work 5] we have used a 
model with two discrete levels (Fig. For discrete 

levels the continuous dynamics, Eq. (^), is replaced by 
jumps between the potential minima, and the appropri- 
ate rate equation in the presence of the external ac field, 
h {t), can be solved by using the adiabatic approximation. 
In such a manner we have accounted for the main exper- 
imental results such as a set of maxima of the output 
signal as a function of H for different fixed ac strength 
ho 0. 

In the experiments analyzed in the previous publi- 
cation the external ac field was the sum of three 
monochromatic waves 



h{t) = hci cos(ajt) -|- 0.45/io[cos(tj 
+ cos(w — il)t] 



(5) 



one with a frequency u, and two with satellite frequencies 
uj ztil, where $7 ^ w. This multifrequency excitation of 
a non-linear system allows one to detect a signal at low 
frequency fl. 

Unlike the previous work, in the new experiments de- 
scribed in the previous sections , an additional pump- 
ing field hp{t) — h2ficos{2flt) was applied to the sys- 
tem. Substituting this field together with Eqs. H2I4|) into 
Eq. ^ one gets 



+ho cos ujt + 0.45/io[cos(ijj -I- Q)t 

+ COs(ljJ — n)t] + /l20 COS 2Ht 



(6) 



Since <C one can decompose AI (t) into a sum of 
two terms as follows 



M (t) = a: (t) + ^ sin (ujt) + 



(7) 



The first term on the right-hand side of Eq. ||7J) will be 
assumed to vary significantly only over times of the order 
of t, while the other terms vary rapidly. Substituting 
into © one can perform an averaging over a single cycle 
time of sin (ujt) . All odd powers of sin (cut) vanish under 
the average while the sin^ (ujt) term will give i. Finally, 
one obtains the following equation for A (t) , the mean 
value of M (t) during the oscillation period, A{t) = < 
M{t) >, 



^ - A[a-bcos (fit) 



(8) 



ccos{2nt)] 
+p'^A^ + /3 = cos {2nt) 

where we have introduced the dimensionless parameters 
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FIG. 6: (Color online) Experimental data redrawn in the di- 
mensionless amplitudes of the output and the input signals, 

4^ and ^^r^ for the constant value of the dimensionless pa- 
no "0 

rameter „ „ . 
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It is clear now that for dimensionless amplitudes of the 
output and the input signal, ^ and !^ , the main param- 
eter which defines the input-output relations is (^-jf^^rff^ 
while the phenomenological parameter /3 is less impor- 
tant. 

Based on the theory described above we may now re- 
draw the experimental data shown in Figs.|31and^ In or- 
der to straighten out these experimental data one has to 
use the scaled dimensionless coordinates and 4?^ , and 
grouping together the experimental points correspond- 
ing to the same value of the parameter ^ ^'^ . These 
regrouping data from Figs. Eland 01 are show in Fig. (HI 

As one can see from this figure the experimental data 
are now well ordered as they should be according to 
Eq. (jSJ. The small differences between different curves 
might be attributable to the different values of parameter 
p. The regrouping procedure was used for the evaluation 
of Hc3 from the experimental data as well. The universal 
behavior Aa(/i2n) occurs possible for different tempera- 
tures. Fig. [Tldemonstrates this type of universal behavior 
for Afi(h2n) in dimensionless coordinates at T = 7.5 K 
and T = 6 K . 

The universal form of the curves in Figs. 16171 suggests 
the validity of our theory. Eq. (jSJ describes the over- 
damped Duffing oscillator with parametric and forced 
periodic excitations. Although an extensive literature is 
devoted to an investigation of the Duffing oscillator (^]), 
the periodic excitations were analyzed only for the un- 
derdamped case (0). We have performed comprehen- 
sive numerical integrations of Eq. 0. A typical result 
of these calculations for some values of the parameter 
^ is shown in Fig. 

Several conclusions can be reached from the analysis 
of this and similar graphs : 

1. Over a wide range of the parameter ^ ^"^h ■ tl^^re 
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FIG. 7; (Color online) (a) An{h2Q dependencies for T — 7.5 
K and T = 6 K for the constant value of the dimensionless 
parameter ^ '^'Lu ■ (b) The same experimental data redrawn 
in the dimensionless amplitudes of the output signal and the 
parametric pumping amplitude, ^ and 



FIG. 8: (Color online) Two bifurcations of the amplitude of 
the output signal An as a function of the amplitude h2n of 
the pumping signal. 
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are two bifurcations at small ("left") and large ("right" 
) values of the amplitudes of the pumping field. 

2. With a decrease of ( g ^^le coefficient which 
multiplies the nonlinear term, the right bifurcation oc- 
curs at a larger amplitude of the pumping signal ,and 
the jump becomes larger. 

3. With a decrease of b (and c) which describe the para- 
metric excitation, the left bifurcation occurs at a some- 
what larger pumping amplitude, and the jump becomes 
slightly larger. 

The currently available experimental data show at 
least qualitative agreement with the suggested theory: 

a) Comparison of the experimental graphs in figures |3| 
and 21 and their scaling treatment shown in Figs. I6I7I 
clearly proves that the theoretical model addresses the 
main features of the experimental data. 

b) The giant increase of the output signal in the pres- 
ence of the pumping field shown in Fig. [21 is present in our 



FIG. 9: (Color online) Experimental data describing the de- 
pendence of the amplitude of the output signal An as a func- 
tion of the amplitude of the pumping signal h2n for different 
values of H. 



theoretical model as well. For instance, for p = 2 * 10~^ 
, equation shows a two order of magnitude increase 
in An in going from = to /120 = 50. 

c) The presence of the left bifurcation can be seen from 
experimental graph which shows as a function of 
h2o for different values of H. The values of /i2n at max- 
ima are shown in Fig. 1101 for different H. In the same 
figure we show also the same quantities obtained from 
Eq. (jSJ. Since in the theoretical analysis we omit some 
constant coefficients (T = l,p « {Hcs — H)~ , etc.), in 
the comparison of theoretical results with experimental 
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FIG. 10: /i2n of maxima obtained from experimental graphs 
Fig. El and from the theoretical analysis of Eq.|Hl 



the right bifuracation, one needs the higher amphtudes 
of the pumping signal which cannot be achieved with the 
existing experimental set up. 



V. SUMMARY 

We have performed an extensive experimental study of 
the non-linear response in surface superconducting states 
of single crystal niobium. Application of a pumping ac 
magnetic field results in a giant (about three order of 
magnitude) increase of the amplitude of the output sig- 
nal. The theoretical analysis, based on the extend double 
well potential model, provides a qualitative explanation 
of the experimental results as well as predicting a new 
phenomenon, the abrupt change (bifurcation) of the non- 
linear response in ac driven superconductors. The latter 
effect, in its turn was partially supported by experiment. 



data one has to set one of the points on the theoreti- 
cal curve with one on the experimental curve, which has 
been done in Fig. ^1 

Both curves are almost straight lines with the correct 
sign of the slope, i.e., we conclude that our theoretical 
model gives at least semi-quantitative agreement with 
experiment. However, one cannot see the right bifur- 
cation on the experimental graph 8 which is, probably, 
connected with some ignored factor in our theory which 
smears this bifurcation. Another possible explanation 
lies in the fact that in order to see in the experiment 
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